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^ , Abstract 

We interpret the dijet enhancement reported by the ALEPH collaboration in the 
process e~^e~ — >4 jets as being due to the production of a pair of bottom squarks, 
r~| ! followed by their R-parity violating decays into pairs of light quarks. Constraints 

on this speculative interpretation are examined. Some of the consequences of our 
, ■ ; ■ hypothesis are drawn within the context of softly broken supersymmetry. 
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Recently the ALEPH collaboration presented a preliminary analysis of about 6 pb^^ of 
data collected at i/i = 130 — 136 GeV at LEP 0]. Although cross sections for standard 
processes appear to be consistent with expectations, they reported some peculiar 4-jet events 
in their data for which there is no canonical interpretation. What ALEPH sees is an excess 
of 4-jet events (14 observed, 7.1 expected), with 8 of these events clustered at a dijet "sum 
mass" of about 110 GeV. This dijet sum mass is arrived at by combining together the 2 dijet 
masses in the events which have the smallest mass difference between them. Assuming that 
what is seen is a signal for something, it corresponds to a cross section of 2.5 ± 1 pb. This 
"signal" is quite distinct from what is expected from QCD 4-jet events, where the event 
distribution is essentially flat in the dijet sum mass. Indeed, in the sum mass bin from 
102-116 GeV where the 8 events are clustered, one expects only 1.35 events from QCD. 

ALEPH makes no claims about this signal and it could well be a statistical fluctuation 
which will disappear as more data is collected. Nevertheless, since the dijet mass difference 
in the analysis is restricted to be below 20 GeV, it is tempting to speculate that the ALEPH 
excess is due to the production of a pair of particles of mass of about 55 GeV which then 
each decay into pairs of jets. In fact, already Farrar ^ has suggested that the ALEPH 
events are associated with the pair production of squarks, which subsequently decay into a 
quark and a light gluino. Although the light gluino scenario is interesting, here we wish to 
speculate in another direction. We also want to associate the ALEPH signal with the pair 
production of squarks. However, in contrast to Farrar, we suggest that what is produced 
is only the bi and that the dijets result from the R-parity violating decay of this squark to 
pairs of light quarks. In what follows we will try to justify this speculation and draw some 
of its consequences. 

For a 55 GeV squark at LEP 1.5 energies, the cross section for e~^e~ —>■ bibi is about 
1 pb. Thus the size of the signal is of the right magnitude. Although the electroweak 
symmetry does not permit an R-parity violating trilinear term involving the {t,b)L doublet 
in the superpotential, there is a possible AB = 1 R-parity violating term involving bpt 0. 
Since the bi squarks mix with the bR squarks, these latter couphngs (if they are present) 
allow for the decay of 6l into light quarks. One must check, however, that this decay does 
not run afoul of any of the existing constraints on the trilinear R-parity violating couplings. 

There are nine possible AB = 1 trilinear R-parity violating couplings in the superpoten- 
tial involving right-handed isosinglet quarks 

W = \,,D^^D^^Ul (1) 

since \ijk = —Xjik- The most stringent bound on these couplings is that on A^^^, which 
comes from the process NN —>■ KKX, where is a nucleon and i^' is a particle with unit 
strangeness. In Ref. @], it was shown that one could have a bound as strong as \ Xdsu\ < 10""^, 
but with large hadronic uncertainties and significant model dependence. Six of the couplings 
Xijk will involve a bR. Xmu\ Xbsu] Xmc] Xbsc] Xbdu Xbst- The last two couplings are not relevant 
here since the top is too heavy to be produced in the decay of a 55 GeV b^. However, 
assuming that 6l is 55 GeV, one must make sure that 

Xbdu Xbst << e 
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to preserve the decay t Wb as the main top decay mode. To our knowledge, there 
are no stringent bounds on Xbsu, \dc and Xbsc- However, \mu gives a contribution to n — n 
oscillations and the experimental limit on the neutron oscillation lifetime of r > 1.2 x 10^ 
sec bounds this parameter to be below about 10^^ [|. As long as some of the couplings \ij 
are of order 10~^ or larger, the decay of the produced bi into light quarks will not produce a 
displaced vertex 0. Thus, the existence of at least one R-parity violating coupling of sufficient 
strength among Xhsm \dc and X^sc makes the suggested scenario phenomenologically viable. 
It could be argued that we have introduced unnatural hierarchies in the couplings Xijk 
by requiring that we evade certain bounds while still retaining one coupling of sufficient 
strength. However, it is worth pointing out that even within the minimal standard model 
one has sizable hierarchies in the Yukawa sector, such as m^j/m^ ~ 10~^. We should note 
that, at this stage, one could have also imagined that the ALEPH events came from the pair 
production of the lightest stop, followed by their R-parity violating decays into light quarks. 
For this to be viable, however, one would have to imagine that the coupling Xsdt dominated 
over both X^st and to account for the lack of decays with a 6-jet in the final state in the 
ALEPH 4 jet sample 0. 

To proceed, however, we must still check whether our suggestion is theoretically sound. 
For the decay 6^ 2-jets to dominate, it is important that the bi be the LSP. In particular, 
all neutralinos should be heavier than the 6/,. Otherwise the decay bi — > bx\ is likely 
to be dominant, since the relevant coupling is of 0(e). We have associated the ALEPH 
enhancement to the production oi a.bi because this squark, along with the stops, has a mass 
that is sensitive to the large top Yukawa coupling. Thus, as is well known 0, even starting 
with universal SUSY-breaking scalar masses at the GUT scale, the bi has the possibility 
of obtaining a nonuniversal mass. In fact, assuming universal soft breaking masses for the 
scalars and the gauginos and tan/5 ~ 2, one finds 



while 



- "^0 + 6.8m?/2, (2) 

m?^ ^ O.Slmo + 5.5mi/2- (3) 

One sees from the above that it is possible to have the 6^ mass be smaller than that of the 
other down squarks, provided the universal scalar mass mo dominates over mi/2- However, 
because we want irti^ ^ 55 GeV then both mo and, particularly, the universal gaugino mass 
mi/2 must be quite small. This, in general, is unacceptable since it leads to one or more very 
light neutralinos in the spectrum. Neutralinos in this mass range are excluded experimentally 
by LEP 10]. But more importantly, the presence of neutralinos lighter than the bi would 
alter the decays of these squarks in an undesirable way. Indeed, in the presence of such light 
neutralinos, the weak decay 6l Xi^ followed by the i?-parity violating decay Xi ~* 3 j^ts 
would give a different experimental signature from that suggested by the ALEPH data. 

^The actual bound depends in detail on the supersymmetric spectrum, as discussed by Goity and Sher 

®- ~ . . . . 

^ In fact, the decay of the via the i?-parity violating operator (|l|) can only occur in the presence of 

some (typically small) bL-hn mixing. Hence we need Aby sin6'b > 10~^, where sin6';, ^ 5 x 10~^ is the sine 
of the bL-bn mixing angle. However, the combination Xhij sin 0b also cannot be too large, for this would lead 
to Z ^ qq'bL decays at an unacceptable level. 
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The presence of light neutrahnos, however, is a feature of the particular pattern of SUSY 
breaking which one assumes at the GUT scale. The principal contributor to the mi/2 piece 
of in Eq. (|^) is the gluino component. However, the neutralino masses are sensitively 
dependent on the SUSY breaking masses one gives to the electroweak gauginos, but 
weakly dependent on the gluino mass. If the electroweak gaugino masses are taken to be 
different from the gluino mass, it is possible to make the neutralinos sufficiently heavy 
by having the SU{2) x f/(l) gaugino masses mi and m2 heavy. This, per se, should not 
affect terribly the mass. However, to be sure one must examine anew the spectrum of 
supersymmetric excitations one gets in the case where the soft breaking of supersymmetry 
involves nonuniversal gaugino masses, mi 7^ m2 7^ 0. 

We present below the results of a study of a minimal supergravity model with unequal 
soft supersymmetry breaking mass parameters in the gaugino sector, but where one still has 
a common SUSY breaking mass for all the scalars. To constrain the model further, we also 
assume that the electroweak symmetry is broken radiatively ||^. Even though a common 
mass is assumed for all the scalars at the GUT scale, these masses evolve to different values 
at low energy as a result of radiative effects. For the squarks of the first two generations 
and for the sleptons, the effective masses at low energy are sensitive functions of the gaugino 
masses and the evolution of the coupling constants. The masses for the third generation 
squarks depend, in addition, on the top Yukawa coupling and its evolution. Furthermore, 
for the mass squared of the and tn one cannot neglect the SUSY-preserving contribution 
of ml, due to the large top mass. 

Solving the renormalization group equations for the squark and slepton masses 0, one 
obtains the following approximate formulas for the light squarks and the sleptons: 

ml - 0.27 cos 2/3M| + xi{t)ml + X2{t)ml 

mr, + 0.5 cos 2/3M| + a;i(t)mi + X2{t)ml 

ml - 0.23 cos 2/?M| + 4xi (t)m2 

ml + 0.35 cos 2/5 M| + -xi{t)m\ + X2{t)ml + x^{t)ml 

9 

ml - 0.42 cos2/?M| + ^Xi{t)ml + X2{t)ml + x^{t)ml 

16 

ml + 0.15 cos 2/5 M| + — Xi(t)m? + X3(t)m^ 

9 
4 

ml - 0.08 cos 2/3M| + -Xi{t)ml + a;3(t)m^ . 

9 

Here t = In and the Xi{t) are functions that depend on the running of the standard model 
coupling constants. Taking Mx = 2 x 10^^ GeV and acuT = jts approximately, 

a;i(t) = 0.038; X2(t) = 0.49; a;3(t) = 6.30. 

Nonuniversal gaugino masses have been advocated recently by Roszkowski and Shifnian|^ in a different 
context. 
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For the third generation squarks, except b^, there are additional contributions due to the 
top Yukawa couphng. One finds: 



m?^(t) = mj^^{t) + hQ{t)ml + hij{t)mimj 
mf^{t) = ml^{t) + 2ho{t)ml + 2hij{t)mimj + 171^ 
m}^{t) = m\^{t) + ho{t)ml + hij{t)mimj + . 

Here the functions /io(^) and hij{t) depend on the evolution of the top Yukawa coupling. For 
example, using as input tan/5 = 2.4 and nit ~ 180 GeV, one has 

ho{t) = -0.494 



hit) 



-0.0154 -0.0004 -0.0024 
-0.0004 -0.1214 -0.0203 
V -0.0024 -0.0203 -1.0793 



Because ho{t) is negative, the scalar mass contribution for Bl, and tn is smaller than 
that for the other squarks. Furthermore, because m| is very weakly dependent on m\ and 

hii < 0, one sees that large values of the f/(l)-gaugino mass will further reduce the 6^ mass 
provided that the contributions of the other two gaugino masses are relatively contained. 
This amounts to a considerably fine-tuned cancelation of the U{1) gaugino contribution 
against the contributions of the other gauginos. A large value for mi has the desired property 
of raising the neutralino masses, securing the role of bi as the LSP. Taking, for example, 
niQ = 20 GeV and 

mi = 1900 GeV; ma = 160 GeV; mg = 80 GeV 
gives mr ~ 56 GeV, which is more than 50 GeV below the next two lightest squarks 



m 



115 GeV and m^^ ~ 260 GeV). [] 



Once all the SUSY-breaking parameters m^ are fixed (along with the SUSY-preserving 
mass parameter /i, which is then a function of tan/?), it is straightforward to deduce the 
gluino mass and the neutralino and chargino mass spectra. For the parameters detailed 
above, one finds: 



m-g = 230 GeV, 



"'it is not necessary to assume a value for tan (3 if one assumes that the top Yukawa couphng at the 
top mass is that which corresponds to the IR fixed point. Ahhough the value for tan/3 we shall use gives 
a Yt{mt) close to the fixed point value, we shall take tan/3 as a free parameter, and fix the top Yukawa 
coupling using the known top mass 180 ± 12 GeV pl . 

^ Here ii is the lighter of the two stops and the detailed value of its mass depends on the mixing between 
and in. This in turns depends on the value of the trilinear SUSY-breaking parameter at the GUT 
scale and on the supersymmetric mass parameter /i. We have taken Aq — and have fixed ^ from the 
minimization condition which must be satisfied to obtain SU(2) x U{1) breaking pO| . For this example one 
has /i = 302 GeV. 



4 



while 



m,.o = 106 GeV; m..o = 305 GeV; m..o = 327 GeV; m.o = 788 GeV 



X2 



X3 



and 



m+ = 106 GeV; m+ = 332 GeV. 



For these same parameters the tree level Higgs masses are 



ruH, = 64 GeV; nin^ = 532 GeV; 



niA 



528 GeV; mH+ = 534 GeV. 



To get an idea of how "fine tuned" the above mass spectrum is, we present in Fig.l a scatter 
plot of minimal neutralino and chargino masses obtained by varying the input parameters 



and tan P, but requiring still that m^^ ~ 55 GeV. 



By lifting the assumption of universal gaugino masses at the GUT scale, we have suc- 
ceeded in producing a spectrum of SUSY particles which exhibits the essential features nec- 
essary to account for the ALEPH "signal" : the LSP is a bottom squark of mass ~ 55 GeV 
which decays to pairs of jets via i?-parity violating couphngs. However, to establish the phe- 
nomenological viability of this scenario, we must still examine various indirect constraints on 



the model. Two such constraints are provided by the p-parameter |12[ and SUSY-mediated 
fiavor changing neutral currents. A further set of constraints can be obtained from analyzing 
new top quark decay modes that are present in this scenario. 

The supersymmetric contributions to the p parameter have been studied previously in 
Refs. O, 111]. It has been shown that the extra Higgs particles and the gauginos of the SUSY 



standard model give a negligible contribution to p. However, non-degenerate SU{2)l squark 
doublets can give a large contribution to p, in complete analogy with the similar result for 
quarks in the minimal standard model. In particular, for the spectrum described above, 
the large splitting between the top and bottom squarks can give a sizable contribution to p: 
writing p = l + SpusM + ^Psusy, we have |l|, |l4l 



Spi 



3a 



SUSY 



SttM^ sin^ Ov. 
cos^ 9t (^cos^ 9b f {mi^ , J + sin^ 9b f {m^^ , m^^ ) 

+ sin^ 9t (^cos^ 9bf{mi^ , "^fej + sin^ 9bf{mi^ , m^J 

- cos^ 9t sin^ 9tf{rni^ , m^-J - cos^ 9b sin^ 9bf{mi^ 



where 9t is the mixing angle between and ti, 9^ is the mixing angle between and 6i, 
and f {1711,1712) is given by 



/(mi, ma) 



171 



22 2 
mimi , mf 

log^. 

mo 



ml 



5 
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lightest chargino mass 



140 



Figure 1: Lightest chargino and neutrahno masses obtained for a sample of 1000 spectra 
with m^^ = 55 ± 5 GeV. 
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For a top mass of 180 ± 12 GeV, a fit to the data yields a liberal f\ bound on Spsusv of 0.004. 
Evaluating SpsvsY for the spectrum given above, we find that the squarks give a significant, 
but not unacceptable, contribution to the p parameter: SpsvsY = 0.0028. 

Supersymmetric particles with the mass spectrum considered here can also have effects 
on flavor changing neutral currents in the B system. In particular, it is known that su- 
persymmetry can give dangerously large enhancements of the 6 — >■ 57 decay rate. The 
branching ratio for charmless radiative B decays has been measured by CLEO |]TB[; at the 
95% confidence level, they have reported 

1 X 10"^ < BR{B X,7) < 4 X 10"^ 

The contribution of supersymmetric particles to this decay rate can be computed using the 
results of Refs. |]T^ |T^, |T^, |2D|. There it was shown that 



|C7(/^)P 



where g^nic/mb) is the phase space factor for the semileptonic decay, and Cj{p) is the 
coefficient of the flavor-changing operator 



BR(5 - 


-^s7) _ 


\v:sVti, 


^ 6aQED 


BR(5 - 




\vl\ 


7rg{mjmb) 



07 



47]-2 



mbSLCTuubRF 



evaluated at a scale p ~ rrib- The QCD evolution of Ct{p) has been computed in Refs. [|T^ [18 



and the contributions of supersymmetric particles to the relevant Wilson coefficients can be 
found in Ref. In the leading logarithmic approximation, we have |T9 | 



C^{p)c^7]-^C^{Mw) + 



14 
Tj 23 



16 
Tj 23 



i=l 



where Cg is the coefficient of the chromo-magnetic moment operator 

9s 



47r2 



and 7] = as{Mw)/<ys{p)- The coefficients hi and Oj are pure numbers independent of any 



model parameters, and can be found in Ref. |]T9 



The coefficients C7{Mw) and Cs{Mw) can be found in Ref. 0]. Normalizing them 
appropriately, the standard model contributions are given by 



and 



3x^ - 2x^ 
4(x- 1)4 



logo; 



5X 



Sx^ + Ix 



-3x2 
4(x- 1)4 



logx + 



24(x- 1)3 
—x^ + + 2x 



(x-1) 



^By this we mean that this is the maximum allowed value of 6p when its correlation with the Peskin- 
Takeuchi S parameter is taken into account, and S is allowed to vary such that the 90% confidence level 
limit on (5/OsuSY is maximized. If we fix S* = 0, the bound decreases to 0.003. 
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where x = mf/M^^. In the present case, we anticipate that the charged Higgs, the charginos, 
and the gluino will give the bulk of the supersymmetric contribution. The charged Higgs 
has long been known to give a sizable enhancement of radiative B decay rates ||T^. The 
charginos can further enhance or suppress the rate [pTI]. The presence of relatively small 



gluino and masses in the spectrum given above indicate that the gluino may also give a 
sizable contribution. The contributions of the charged Higgs and the charginos to Ci^s{M\y) 
may be found, for instance, in Ref. |2^. For the gluino, we have pi 



C|(M, 



a. sm 



'9 



w I 



a 



2 

bL 



4^2 



9(z 



log z + 



-4^2 _ + 2 



27(^-1) 



and 



C|(M, 



a, sm' 



w 



sbL 



2a 



VtbVtl 



■9z + , 



-11^2 + 40^ + 19 



18U- 1 



where z = m'^/m'^ , and Ugg'^ is the flavor mixing matrix that appears in the coupling of a 
quark to a gluino and a left-handed squark. We have neglected a small additional contri- 
bution that involves mixing. Similar "super-CKM" matrices appear in the couplings 
of the charginos to bottom quarks and top squarks. Given the large non-degeneracy of the 
squarks in the above spectrum, it would be reasonable to expect that these "super CKM" 
matrices may deviate from the standard model CKM matrix. Hence in the following we will 
retain these CKM elements to make explicit the dependence on these angles. 

Evaluating Cy^s for the spectrum given above (m^^ = 56 GeV and rrig = 230 GeV) and 
including the contribution of the charged Higgs and the charginos, we find an estimate for 
r(6^ S7): 

r(6 ^ S7)[SUSY1 , 

-^^ ^{[^ J ~ 1 + 0.18 + 0.09a - 0.49A 2 

r(6 S7)[MSM] ' ' 

where the four terms come from W, charged Higgs, gluino, and chargino loops. The contri- 
butions of all of the squarks have been included and expressed in terms of third generation 
mixing angles using CKM unitarity. The quantities a and A are given by 



a 



U- TI*- 

^bbL^sbL 

VtbVtl 



and 



A 



VtbV,l 



The matrix V is the super CKM matrix describing couplings of charginos to up-type squarks 
and down-type quarks. We see that for A ~ 1 (a reasonable value), the charginos interfere 
destructively with the W and charged Higgs, reducing the rate. The possibility of such a 
phenomenon has been noted in Ref. [E^. Given the standard model estimate pOl 



BR(5 ^ X,7)[MSM] ~ 1.9 ± 0.5 x 10' 



we conclude that the rate is compatible with the CLEO determination. 
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The presence of light supersymmetric particles in the spectrum allows the top quark to 
decay in other modes besides the standard t Wb mode. For the model spectrum discussed 
above, the dominant non standard decay is t — > bixt Q- The branching ratio for this mode is 
of 0(30%) which, although sizable, is probably acceptable given the experimental uncertainty 
in the top cross section |TT|. Because the dominant chargino decay is xt ^2^, with the 
bi decaying into dijets, the final state for these non standard top decays will contain five 
jets. This is not a particularly easy signal to dig out. Nevertheless, if some top tagging can 
eventually be implemented at the Tevatron, looking for multijet decays of the accompanying 
t may be the best way to dig out the b^ in hadronic interactions 0. 

It is quite likely that, in the end, the ALEPH dijet enhancement will prove to be just 
a statistical fluctuation, making the scenario discussed here moot. Even if this were to 
turn out to be the case, we believe that some of our disquisitions may continue to prove 
useful. Three lessons which emerge from our speculations stand out in particular. First, 
significant deviations from the spectrum predicted by minimal versions |T^ of softly broken 
supersymmetry can occur as a result of some simple changes in the underlying assumptions 
(e.g. having non-universal gaugino masses). Second, a low mass sbottom (or a low mass 
stop) should not be unexpected, given the large top Yukawa coupling. Third, although there 
is a natural prejudice against R-parity violating couplings, their presence at some level is 
by no means ruled out. If present, such couplings totally alter the "standard" signals of 
supersymmetry. 
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^ Another allowed decay is t ^ tibi^b, but this decay is strongly kinematically suppressed. There may 
also be R-parity violating top decay modes into bd and 6s, but the relevant branching fractions for these 
modes depend on the unknown couplings Xut- 

^Although bL — b^ pairs are copiously produced at the Tevatron, because of the large QCD background 
it appears very difficult to detect them with an ALEPH-like analysis. 
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